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ABSTRACT 

Small-scale magnetic field concentrations (magnetic elements) in the quiet Sun are believed to contribute to the energy budget of the 
upper layers of the Sun’s atmosphere, as they are observed to support a large number of MHD modes. In recent years, kink waves 
in magnetic elements were observed at different heights in the solar atmosphere, from the photosphere to the corona. However, the 
propagation of these waves has not been fully evaluated. Our aim is to investigate the propagation of kink waves in small magnetic 
elements in the solar atmosphere. We analysed spectropolarimetric data of high-quality and long duration of a photospheric quiet Sun 
region observed near the disk center with the spectropolarimeter CRISP at the Swedish Solar Telescope (SST), and complemented by 
simultaneous and co-spatial broad-band chromospheric observations of the same region. Our findings reveal a clear upward propaga¬ 
tion of kink waves with frequency above 2.6 mHz. Moreover, the signature of a non-linear propagation process is also observed. 

By comparing photospheric to chromospheric power spectra, no signature of an energy dissipation is found at least at the atmospheric 
heights at which the data analysed originate. This implies that most of the energy carried by the kink waves (within the frequency 
range under study < 17 mHz) flows to upper layers in the Sun’s atmosphere. 
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1. Introduction 


Magneto-hydrodynamic (MHD) waves are ubiquitous in the so¬ 
lar atmosphere. Both large-scale magnetic structures, extending 
over several Mm in length, and small-scale magnetic concentra¬ 
tions (magnetic elements) at scales close to the resolution limit 
of modern solar telescopes (~ 120 - 150 km on the solar photo- 
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immediately recognized as fundamental for the energy budget 
of the solar upper atmosphere. This hypothesis became even 
more attractive after the availability of very high-resolution 
spectropolarimetric data that have provided strong observational 
proofs that magnetic ele ments also cover a significant fra ction 
of the solar photosphere (lLagg et al.ll2010l: iBonet et al.ll2012l) . 
The ub iquitous although non-homogeneously distributed 
dst angalinii 12014 )) magnetic elements are perturbed by the 
external photosphe ric plasma flows and are advected in the 
solar photosphere (lAbram enkoet all 120111; iLepreti et al.ll2012t 
iGiannattasio et al.ll2013l 2014allbl) . Recent observations of the 
photosphere have revealed that different types of waves can co- 
exist and interact in the s ame magnetic element (Stangali ni et al.1 
201 3bj). iJess et all (120121) reported of propagating compressible 
waves in magnetic elements with periods of 110 - 600 s in both 
the ph otosphere and the chromosphere. Moreover, iMorto netalJ 
(201 3) have observed torsional modes (lErdelvi & Fedunll2007l) 
generated by vortices in the solar photosphere. These are just 
two examples demonstrating the large variety of waves that have 
been observed so far in magnetic elements. 


[Morton et al.ll ( 2014 ) have investigated the propagation of 
incompressible waves to the chromosphere, finding a good 
agreement between photospheric and chromospheric velocity 
power spectra at frequencies below 8 mHz. 

Among the many MHD waves that magnetic elements can 
support, kink waves induced by the solar convection have been 
longly considered as a viable mechanism for transferring energy 
from the photosphere to the upper layers of the Sun’s atmosphere 
jKalkofen 1997; Hasan et al. 12003 : [Musielak & Ulmschneiderl 
200 3bllal : lHasan & van B allegooiien 20081) . 

Kink waves have been reported at different layers in the solar 
atmosphere; from the lower photosphere (e.g. Keys et ahll201 ll) . 
to th e chromosphere (e.g. Jafarzadeh et al.ll20 1 3h and the corona 
(e.g.lMcIntosh et akll201 1 ). 

[Stangalini et al.l ( 20141) have observed the presence of subhar¬ 
monics, with a fundamental period of 7.6 min, consistent with 
the granular time scale. This constitutes an observational evi¬ 
dence for the excitation of kink waves by the solar convection. 
The same authors also argued that the detected subharmonics 
might be the signature of chaotic excitation. Indeed, subhar¬ 
monics represent a general characteristic of complex chaotic 
systems, in the presence of strong forcing (Sander & Yorkel 
l2009ll20T0l) . 

By using both simulations and observations, iKevs et aTI (1201 ll) 
found that the majority of magnetic elements have horizontal 
velocities between 0 and 1 km/s, and 6% of them have veloc¬ 
ities in excess of 2 km/s. The same authors argued that this 
significant fraction of magnetic elements may contribute to heat 
the upper atmospheric layers, since velocities over 2 km/s can 
effici ently drive kink waves (IChoudhuri et al.l Il993). More on 
this, [jafarzadeh et ak| (1201 3l) have observed very short duration 
velocity pulses up to 15 km/s in the internetwork chromospheric 
magnetic elements that are well in excess of their mean value 
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Fig. 1 . Example of the data analysed in this study. Left: photospheric circular polarization map saturated between -3% and 3%. 
Center: chromospheric broad band image of the same region. The 35 longest-lived magnetic elements used in the analysis are 
highlighted by yellow circles. The red crosses mark the six elements whose power spectra are shown in Fig. 0 (see text for more 
details). Right: continuum image taken at 395.3 nm with 1A bandpass. 


(~ 2 km/s), and that may significantly contribute to the energy 
budget of the upper layers of the Sun. 

In addition to the aforementioned detection of subharmonics, 
IStangalini et al.l ( 2013ah observed the presence of high fre¬ 
quency peaks of up to 10 - 12 mHz in the power spectra of 
kink oscillations in the solar photosphere. These high frequency 
oscillations are well above the cutoff frequency expected for 
kink waves, which is always smaller than the acoustic cutoff 
at ~ 5.3 mHz (Spruit! 119811) . Therefore, they are expected to 
propagate upwards. Although the physical properties of the 
magnetic elements analyzed were quite similar, IStangalini et alj 
d2013al) found that each individual magnetic element had its 
own signature in the power spectrum. This demonstrates that 
the kink oscillations are largely influenced by the ambient 
surrounding conditions. 

iMorton et all ( 2014 ) have also studied the average power spectra 
of kink waves in both the chromosphere and the corona. Their 
main finding is that a considerable part of energy is dissipated 
between the chromosphere and the corona itself. This was 
proved by the observed damping of the coronal power spectra, 
although the chromospheric and coronal data were not sampled 
simultaneously, and were not co-spatial. 

However, this study, as well as those reporting the presence of 
kink oscillations at different heights in the Sun’s atmosphere, 
suffers lacking of information on the propagation of these 
waves between different atmospheric layers, due to lack of 
co-spatial and simultaneous data from different heights in the 
solar atmosphere. 


To overcome this limitation, we analysed high-cadence and 
high-resolution simultaneous data of the solar photosphere 
and chromosphere, and investigated the propagation of kink 
waves in small magnetic elements in the quiet Sun, thus 
complementing previous results reported in the literature. 


2. Data set and methods 


The data set analysed in this work was acquired on 2011 Augus t 
6 with the imager s pectropolarimeter CRISP (IS charm er 2006; 

at the Swedish Solar Telescope (SST, 
and was aimed to the study of a quiet 
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Sun region observed at disk center. It consists of a full-Stokes 
spectropolarimetric time series of spectral scans acquired at 30 
wavelengths in the photospheric Fe I doublet at 630 nm, with a 
spectral sampling of ~0.0044 nm. In addition, simultaneous and 
co-spatial chromospheric broadband images taken at the core 
of the Ca II H line at 396.9 nm were acquired. For the spec¬ 
tral lines the following core-formation heights are assumed: i) 
230 km for the Fe I line at 630.1 nm, ii) 200 km for the Fe I 
line at 630.2 nm dFaurobert et al. 20121), and in) 450 km for the 
Ca II H fine at 396.9 nm ( Jafarzadeh et al.ll2013l) . The observa¬ 
tions started at 07:57:39 UT and lasted for 47 minutes with a 
cadence of the spectral scans of 28 s (100 spectral scans). The 
pixel scale was ~ 0”.059/pixel for the Fe I data and 0”.034/pixel 
for the Ca II H data. The Ca II H and Fe I data were co-aligned 
and remapped to the same FoV of area ~ 53” x 53”. The the¬ 
oretical full-width-half-maximum (FWHM) of the point spread 
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(c) Power spectra magnetic element n.3 
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(g) Power spectra magnetic element n.4 
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(k) Confidence level magnetic element n.5 



Fig. 3. Photospheric and chromospheric power spectra of horizontal velocity oscillations, and confidence level plots obtained from 
a randomization test of the velocity time series. 


function (PSF) at 630 nm is 0”.16 or 120 km on the solar pho¬ 
tosphere. This constituted the diffraction limit of the telescope. 
The spectral resolution of CRISP at these wavelengths is 0.0055 
nm (ide la Cruz Rodriguez et al]l2015b . 

The standard CRISP calibration procedure was applied to 
the dafiL This includes the CRISPRED calibration pipeline 
(ide la Cruz Rodriguez et al.1120151) , and the MOMFBD (Multi 


Object Multi Frame Blind Decomposition, Ivan Noort et al.l 
120051 ) restoration aimed at limiting seeing-induced distortions in 
the images. The photospheric Fe I data were used to estimate the 
total circular polarization CP within each pixel, which is defined 
as 

r dA 

CP = I \V(A)\—, (1) 
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Fig. 2. Chromospheric and photospheric velocity amplitude of 
the 35 analysed magnetic elements, as a function of their mean 
velocity. 


where V(A) is the V-Stokes profile, and I c the average continuum 
intensity. 

We remark that CP signals can be attributed to the line-of-sight 
(LOS) magnetic field. As a consequence, in what follows small- 
scale CP signals can be identified with magnetic elements. 

In Fig. Q] we show an example of signed CP polarization map, 
where the sign is assumed to be that of the blue lobe in Stokes-V 
profiles, and the corresponding simultaneous broadband image 
in the core of the Ca IIH line. In the same image we also show a 
photospheric continuum image of the same field taken at 395.3 
nm (1 nm width). 

Our analysis is based upon the identification and tracking of 
the longest-lived magnetic elements observed at both the pho¬ 
tospheric and chromospheric heights. 

The studied magnetic elements wer e identified and tracked bv 
using the well tested YAFTA code (IWelsch & Lon ticonell2003t 
IDeForest et al.l 120071) . providing it with the CP maps in input. 
The parameters of the code were set to track only those mag¬ 
netic elements i) whose area was larger than 4x4 pixel 2 (or 
equivalently 0”.24 x 0”.24), which takes into account seeing- 
induced adaptive optics (AO) residual aberrations that make the 
FWHM larger than the theoretical (seeing-free) value; ii) whose 
signal in the CP maps was above two standard deviations of the 
CP image. A total of 185 magnetic elements was found to match 
these criteria. Among them we selected the longest lived (at least 
90 steps or equivalently ~ 42 min), and whose chromospheric 
counterpart was clearly identifiable by visual inspection of the 
simultaneous data. 

This yielded 35 magnetic elements, which are marked with yel¬ 
low circles in Fig. |T| 

The mean effective radius of the studied magnetic elements, de¬ 
fined as V Area/n, is ( R e ) - 280 km. 

The magnetic elements selected are, on average, co-aligned 
within ~ 2 pixels. Given the height difference of the sampled 
atmospheric layers (~ 220 km), this yields an inclination an¬ 
gle of ~ 22 de grees. This value is in good agreement with 
iJafarzadeh et al.l?2014l) . The tracking of the chromospheric fea¬ 
tures was performed by applying a procedure based on the search 
for the center of mass of the intensity distribution in windows of 
area 10x10 pixel 2 encompassing the intensity enhancement co- 
spatial to the photospheric feature. 

At each time step, the horizontal velocity of each identified mag¬ 
netic element was computed as the time derivative of the mea¬ 
sured position of the feature. 


We applied a standard FFT-based analysis to compute the hor¬ 
izontal velocity power spectrum of each photospheric magnetic 
element and its chromospheric counterpart. In doing this, only 
one component of the horizontal velocity was considered. This 
is to avoid frequency doublings due to the square root operator 
used in the composition of the two components of the horizontal 
velocity. 


3. Results 

In Fig. [2]we plot the photospheric (black plus symbols) and chro- 
mopsheric (cyan diamonds) RMS (root mean square) velocity 
amplitude of the magnetic elements as a function of their mean 
velocity i.e. the mean velocity during the whole lifetime of the 
element. As we can see, the velocity amplitudes are appreciably 
larger in the chromosphere than in the photosphere. This can be 
also seen in Fig. [3] where we show the power spectra of the hori¬ 
zontal velocity for six representative cases among the 35 studied. 
The chromospheric power spectra are largely amplified with re¬ 
spect to their photospheric counterpart. On top of that, they are 
also characterized by the presence of different peaks from 2-3 
mHz, up to 13 - 15 mHz. It is worth noting that the amplifica¬ 
tion is not constant throughout the whole spectral range. Indeed, 
different spectral bands are more amplified than others. The dif¬ 
ferent magnetic elements do not share the same peaks, and each 
of them has its own spectral signature. 

In Fig. [3 we also show the confidence levels of each pair of 
spectra obtained from a randomization test. 

This test consists in the random shuffling of the measurements 
of the time series associated to the horizontal velocity of each 
magnetic element. More in detail, this shuffling is repeated a 
thousand times for each time series and, at each step, the power 
spectral density is computed. The number of times a certain peak 
in the spectrum is larger than that obtained from the original time 
series (i.e. the correctly ordered one) is counted. This gives the 
probability that a certain peak in the power spectrum is obtained 
from a noisy time series, and constitutes the confidence level. 

As clearly shown by the results obtained, the confidence level 
associated to most of the prominent peaks in the power spectra 
have a very high confidence level (above 0.9). For this reason the 
peaks shown in the power spectra of Fig. [3 can be regarded as 
highly trusted. 

We also studied the wave propagation by applying a phase lag 
analysis. For each magnetic element the FFT phase and co¬ 
herence spectra between the photospheric and chromospheric 
heights were estimated. Only those phase measurements in the 
spectrum whose coherence level is larger than 0.9 are considered 
for the estimation of the global phase spectrum obtained from 
the whole sample. This is necessary to avoid the inclusion of un¬ 
reliable phase measurements in the final global phase spectrum. 
We remark that a 0.9 coherence level represents a very restric¬ 
tive threshold, thus ensuring the reliability of our global phase 
spectrum estimate. 

The results of this analysis are shown in Fig. QJ where we plot 
the unwrapped global phase spectrum obtained by collecting the 
most reliable phase measurements (those with coherence larger 
than 0.9) in the phase spectra among the 35 selected magnetic 
elements. 

In the same figure we also plot for comparison the phase 
spectrum expected for propagating kink waves in thin flux 
tubes as obtained from the dispersion relation of Eq. 4 in 
IPloner. S. R. O. & Solanki, S. K.1 dl997h . The phase relation 
is calculated for a propagation speed of 6 km/s, a pres¬ 
sure height scale of ~ 130 km, and a kink cutoff fre- 
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Fig. 4. Phase diagram obtained from the tracking of the selected 
magnetic elements (diamonds). Among the phase spectra used 
and belonging to the selected elements, only the measurements 
with a coherence above 0.9 are retained. The dashed line repre¬ 
sents the expected phase for propagating kink waves at 6 km/s 
and with a cutoff frequency of 2.6 mHz, between two layers with 
a separation of 220 km. 


quency of 2.6 mHz. These are in fact typical values ex¬ 
pected for small scale flux tubes (see for example ISprui ll98lt 
IPloner, S. R. O. & Solanki, S. K.I11997]) . In our sign convention 
a positive phase difference represents upward propagation. 

We then studied the amplification of the power at chromospheric 
heights. 

The longest-lived 35 magnetic elements were considered to es¬ 
timate the average power spectrum in both the photosphere and 
the chromosphere. In Fig.[5]we show the average chromospheric 
and photospheric spectra (upper panel) and the amplification 
spectrum (lower panel). The amplitude of kink waves is largely 
amplified in the 2-8 mHz band, with a maximum around 4.2 
mHz where the amplification at the chromospheric heights sam¬ 
pled by the core of Ca II H spectral line reaches a factor of five. 
The error bars represent the standard error of the mean. 


4. Discussions and conclusions 

In this work we have studied the propagation of kink waves in 
magnetic elements between the photosphere and the chromo¬ 
sphere. 

The power spectra of the horizontal velocity of different mag¬ 
netic elements in the photosphere show different patterns. This 
is not surprising since the kink mode is largely influenced by 
the ambient plasma surrounding the magnetic elements, which 
affects their oscillatory behaviour. This result confirms previous 
findings bv IStangalini et all (12013a ). obtained from the analysis 
of Hinode data at lower spatial and temporal resolution. 
However, the results presented here also reveal, for the first time, 
that the chromospheric spectra show remarkable differences with 
respect to their photospheric counterparts. The same magnetic 
element exhibits a completely different pattern of peaks in the 
photosphere, with respect to the chromospheric power spectrum. 
This is an evidence for a non-linear wave propagation regime. 

It is worthwhile noting here that a randomization test demon¬ 
strated the reliability of the obtained power spectra with a confi¬ 
dence level larger than 90%. We remark that the randomization 
test also validates the quality of the tracking. More in detail, it 
is unreasonable to find high confidence levels in the spectra of 


noisy signals. For this reason, the high confidence levels found 
in the spectra, together with the large coherence found between 
the two sampled layers, can be regarded as a good validation of 
the methods employed for the tracking of the magnetic elements 
at both heights. 

We also studied the propagation of the kink waves between the 
photospheric layer sampled by the Fe I 630.1 nm spectral line 
and the chromospheric height sampled by the Ca II H spectral 
data. For this purpose, only the 35 longest-lived magnetic ele¬ 
ments were considered in order to maximize the frequency reso¬ 
lution in the power spectra. 

Moreover, by considering only frequencies in the phase spectra 
with very high confidence level between the photosphere and the 
chromosphere, we estimated the global phase spectrum. This is 
done to ensure the reliability of the phase estimates themselves. 
However, we note here that, since the propagation of the kink 
waves is non-linear, in doing this we are only selecting those 
cases for which a linear regime can be applied. 

The phase diagram so obtained shows the upward propagation 
of the kink waves compatible with a velocity of the order of 6 
km/s, and a cutoff frequency of 2.6 mHz. Such value for the 
cutoff frequ ency is consistent with the theoretical expectations 
(Hi ruit)ll98ll) . and represents a first experimental test of this in 
the lower solar atmosphere. 

It is worth remarking here that the inclination of the magnetic 
elements can significantly affect the theoretical phase relation 
shown for comparison in Fig. [4] This is be caus e the dispersion 
relation of kink waves (see IPloner. S. R. O. & Solanki. S. K.l 
11997! ) depends on the acoustic cutoff frequency, which is in 
turn lowered in inclined magnetic fields (see for example 
iMcIntosh & Je fferies 2006). In our case, the selected magnetic 
elements have an inclination angle of ~ 22 degrees thus they can 
be considered nearly vertical. However, the upward propagation 
of kink waves is unambiguously revealed by the phase analysis 
of Fig. 0 This result is valid independently of the inclination an¬ 
gle and other physical conditions of the flux tubes that can affect, 
to some extent, the exact values of the cutoff frequency and the 
propagation speed (IPloner. S. R. O. & Solanki. S. K.ll 19971) . Our 
findings show that the chromospheric power spectra are largely 
amplified. The amplification spectrum shows that most of the en¬ 
hancement is concentrated in the 1-8 mHz band, although, there 
is no significant power reduction at least up to the limit imposed 
by the data sampling (~ 17 mHz). 

In the presence of a non-linear propagation regime, kink waves 
are expected to dep osit a significant fraction of energy through 
dissipation (iMorton et al.ll2014l) . The amplification spectrum es¬ 
timated from the 35 studied magnetic elements demonstrates that 
this is not the case at the heights sampled by the data analysed 
in this work (up to ~ 450 km above the photosphere). A dissi¬ 
pation of kink waves would in fact result in a deficit of power 
at chromospheric heights, with respect to photospheric ones. For 
this reason, we argue that, in absence of deficits of power in the 
chromospheric spectra, no dissipation should be expected up to 
approximatively 450 km above the photosphere, where the core 
of the Ca_II H spect ral line forms. 

IMorton et al.l(l2014l) have observed the dissipation of kink waves 
between chromosphere and the corona. In this context, the re¬ 
sults presented here demonstrate that the energy contained in the 
MHD kink waves in small magnetic elements is not dissipated 
in the low atmosphere but flows entirely to heights which are 
above the altitude of formation of the core of the Ca II H spec¬ 
tral line. This may have implications for the energy balance of 
both chromospheric and, in particular, coronal heights, to where 
most of the energy flows. Finally, we note that the lack of damp- 
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Fig. 5. Top: Kink photospheric and chromospheric power spectra 
obtained from 35 magnetic elements and averaged in bins of 0.8 
mHz. Bottom: Amplification spectrum. 

ing is also evidence for the fact that the kink waves do not inter¬ 
act with other modes during the propagation between the heights 
sampled by our observations. 
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